Isospin effects on squeeze-out flow in heavy-ion collisions 
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The squeeze-out flow in reactions of ^^''Sn -I- ^^''Sn and ^^''Ba + ^^*Ba at different 
incident energies for different impact parameters is investigated by means of an isospin- 
dependent quantum molecular dynamics model. For the first time, it is found that the 
more neutron-rich system (^^'*Sn + ^^''Sn) exhibits weaker squeeze-out flow. This 
isospin dependence of the squeeze-out flow is shown to mainly result from the isospin 
dependence of nucleon-nucleon cross section and the symmetry energy. 

PACS. 25.70.-z Low and intermediate energy heavy-ion collisions— 25.75. Ld Collective 
flow 

The recent advance in radioactive nuclear beam (RNB) physics provides people a unique opportunity 
to investigate isospin effects in heavy- ion collisions (HIC's) 0J^. Isospin effects on directed flow ^-|^, 
radial flow and rotational flow in HIC's at intermediate energies have been explored theoretically 
and/or experimentally, and it is indicated that the isospin dependence of nuclear collective flow may 
provide people some important information about the isospin-dependent nuclear EOS, particularly the 
symmetry energy, and the isospin-dependent in-medium N-N cross section. The squeeze-out flow |8| is 
of special interest since it comes from the only direction where the nuclear matter might escape without 
being hindered by the presence of the cold spectator remnants and thus a less disturbed information 
on the matter of high density and temperature is expected. Therefore, it is very signiflcant to explore 
isospin effects on the squeeze-out flow. 

In this Short note we report results of the flrst theoretical study on the squeeze-out flow from reac- 
tions of ^^^Sn + ^^"^Sn and ^^''Ba -I- ^^"'Ba at different energies for different impact parameters within 
the framework of an isospin-dependent quantum molecular dynamics (IDQMD) model which includes 
the symmetry energy, Coulomb interaction, isospin-dependent experimental N-N cross sections, and 
particularly the isospin-dependent Pauli blocking ||^,^. In the initialization process of the IDQMD 
model, the neutron and proton are distinguished from each other and meanwhile the nonphysical 
rotations in the initialized nuclei have been removed . In the IDQMD model, the nuclear mean field 
can be parameterized by 

U{p, T,) = a{p/po) + Pip/poV + ^(1 - r^^Wc 

+C^^^^T, + C/^"^ (1) 
Po 

with pq the normal nuclear matter density (here is 0.16 fm^'^); p, p„, and pp are the total, neutron, 
and proton interaction densities, respectively; Tz is the zth component of the isospin degree of free- 
dom, which equals 1 or —1 for neutrons or protons, respectively; C is the symmetry strength; Vc is 
the Coulomb potential; and [/^"^^ is the finite range Yukawa (surface) potential which will vanish for 
infinite nuclear matter. The forms and parameters of Eq. dl|) can be found in Ref. The IDQMD is 
different from the so-called IQMD (Isospin-QMD) model |10|-12| by the Pauli blocking, the initializa- 



tion process, and construction of fragment. This model has been used recently to explain successfully 
several phenomena in HIC's at intermediate energies, which depend on the isospin of the reaction 
system |5|-|7|,p|,|l3| . In the present calculations, the so-called soft EOS with an incomprcssibility of 
K — 200 MeV is used and the symmetry strength C ~ 32 MeV without particular consideration. 

In the QMD model fl^ , the reaction plane is known a priori and it is defined as the x-z plane 
(z-axis corresponds to the beam direction). The azimuthal angle with respect to the reaction plane 
can be written as 
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(j) = arctan{Py/P^), (2) 

where Px and Py are the x- and y-components of nucleon momentum in the center-of-mass (cm.) 
system. Many studies have shown that the azimuthal distribution in HIC's can be fitted well by the 
Legendre polynomial up to the second order jl^ , i.e., 

dN/d(p = c(l + ai cos(0) + 02 cos(2(/))). (3) 

The coefficient ai represents the strength of the in-plane directed flow while a positive 02 reflects the 
strength of the rotational collective motion (the azimuthal distribution peaks at </) = 0° and ±180° 
simultaneously) and negative one (the azimuthal distribution peaks at = ±90° ) the out-of-plane 
squeeze-out. For the mid-rapidity azimuthal distribution, the ratio Rn, 

_ dN/d(j){9Q°) + dN/d(j){-m°) _ 1 - 02 
^~ dN/d,p{0°)+dN/d,p{180°) ^1 + 02' ^' 

measures the strength of the squeeze-out flow in a quantitative way. A value Rn > 1 corresponds to a 
preferential out-of-plane emission. It has been shown that Rn depends on the transverse momenta of 
nucleons and in the present study the normalized nucleon transverse momentum Pt/ Pproj (where 
Pt and Pproj represent the transverse momenta of nucleon and the projectile momentum in the cm. 
system, respectively) is limited to Pt/ Pproj > 0.5 for simplicity. Meanwhile, the mid-rapidity is defined 
as a narrow region around the cm. rapidity by applying the condition —0.25 < (j//2/proj)c.m. < 0.25, 
where (y/yproj)c.m. is the reduced cm. rapidity. In addition, it has been shown that Rn depends on 
the fragment mass p6| and in this work the calculated results are obtained from all nucleons entering 
the analysis to accumulate the numerical statistics. 

Fig. 1 displays the IDQMD model prediction of Rn as a function of incident energy for reactions of 
^^''Sn + ^^^Sn and ^^''Ba ± ^^^Ba at impact parameter h~A fm. The errors shown are the statistical 
errors resulting from the Legendre polynomial fits. It is indicated in Fig. 1 that the Rn value increases 
with increment of the incident energy and then saturate or decrease at higher incident energies, which 
is in agreement with the recent experimental results ||r^ , [l^ . The reduction of the Rn value at 
higher incident energies is easy to understand by the shadowing effects of the cold nuclear matter 
surrounding the participant zone which decrease because the projectile spectators escape faster. A 
more interesting feature in Fig. 1 is that the more neutron-rich system ^^^Sn ± -'^^^Sn displays 
systematically smaller Rn values than the system ^^''Ba ± ^^''Ba, which implies that there exists 
a strong isospin dependence of the squeeze-out flow, namely, the more neutron-rich system displays 
weaker squeeze-out flow. Meanwhile, one can find that the isospin dependence of the squeeze-out flow 
decreases as the incident energy increases, which may be due to the fact that the isospin dependence 
of the JV-JV cross section disappears at higher incident energies. 

In order to investigate further the isospin dependence of the squeeze-out flow, we show in Fig. 2 the 
ratio Rn as a function of impact parameter for reactions of ^^''Sn + ^^^Sn and ^^^Ba ± ^^^Ba at 350 
MeV/nucleon. The errors shown are the statistical errors resulting from the Legendre polynomial fits. 
Similarly, a strong isospin dependence of the squeeze-out flow is observed once again, namely, the more 
neutron-rich system ^■^'^Sn ± ^^'^Sn displays systematically smaller Rn values than the system ^^''Ba 
+ ^^^Ba at different impact parameters. In particular, this isospin dependence is more pronounced in 
semi-central collisions. In addition, it is indicated that the largest Rn value is obtained at 6=6 fm, i.e., 
in semi-central collisions, which is in good agreement with the recent experiment where a maximum of 



Rn located at about 6=6 fm has been evidenced in reaction of Au+Au at 400 MeV/nucleon 17 lq|. 
This phenomenon can be explained by an expansion shadowing scenario where the expansion of 
the participant matter is rescattered by the cold target or projectile spectator. 

It is important to investigate the influence of the symmetry energy and isospin-dependent JV-JV cross 
section on the squeeze-out flow since the isospin dependence of the squeeze-out flow may result from 
the competition among several mechanisms in the isospin-dependent reaction dynamics, such as the 
symmetry energy, isospin-dependent JV-JV cross sections, and so on. 

Using different symmetry energy strength C and parametrizations of JV-JV cross sections, we show 
in Fig. 3 the IDQMD model predicted normalized azimuthal distribution from ^^''Sn ± ^^*Sn (solid 
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circles) and ^^''Ba + ^^"^Ba (open circles) at 350 MeV/nucleon and 6 = 6 fm for mid-rapidity nucleons. 
Meanwhile, the results of Legendre polynomial fits according to Eq. (||) for ^^"'Sn + ^^"'Sn (solid line) 
and -'^^'^Ba + ^^^Ba (dashed line) as well as the resulting oi and 02 are also included in Fig. 3. For 
the results shown in Fig. 3 (a) we use C — 32 MeV and experimental N-N cross section Coxp which 
is isospin dependent. In Fig. 3 (b) we use C = (no symmetry energy) and (icxp- The case of using 
C — 32 MeV and Cugnon's N-N cross section acug which is isospin independent, is plotted in Fig. 3 
(c). 

One can see from Fig. 3 that the azimuthal distribution exists minima at = 0° and ±180° (i.e., 
in the in-plane P^j-direction) and maxima at = ±90° (i.e., in out-of-plane P^-direction which is 
perpendicular to the reaction plane). These features imply that at mid-rapidity more nucleons are 
squeezed out perpendicular to the reaction plane than in the reaction plane. In order to see more 
clearly the isospin effects on the squeeze-out flow, we give in Table 1 the Rn values extracted from 02 
in Fig. 3 for different cases. The errors shown are the statistical errors resulting from the Legendre 
polynomial fits. One can see from Table 1 that both the symmetry energy and the isospin-dependent 
N-N cross section enhance the strength of the squeeze-out flow but the latter enhances it more strongly. 
Particularly, it is indicated that the influence of (Toxp on system ^^^Ba + ^^'*Ba is stronger than that 
on the system ^^''Sn + ^^"'Sn, which is easy to understand since the neutron-proton cross section is 
about three times larger than the neutron-neutron or proton-proton cross section for (Toxp at energy of 
about 350 MeV/nucleon, which results in more iV-JV collisions for ^^^Ba + ^^''Ba. On the other hand, 
the symmetry energy is generally repulsive and the change of symmetry strength C might modify 
the equation of state and thus the squeeze-out flow. From above analysis, one can conclude that the 
isospin dependence of the squeeze-out flow seems to mainly result from the isospin dependence of 
JV-JV cross section and the symmetry potential has less influence on it. In above calculations, only 
the free-space N-N cross sections are adopted and the in-medium effect is only simulated by the Pauli 
blocking. However, the in-medium JV-JV cross sections and their isospin dependence might be strongly 
density dependent ||l9| , po|] . The isospin dependence of the squeeze-out flow may provide a unique 
opportunity to study the isospin dependent in-medium JV-JV cross sections. 

In summary, by using the IDQMD model, we studied for the flrst time the out-of-plane squeeze-out 
flow in reactions of ^^"'Sn + ^^^Sn and ^^^Ba + ^^"'Ba. A strong isospin dependence of squeeze-out 
flow has been found, namely, the more neutron-rich system exhibits weaker squeeze-out flow, which 
is shown to mainly result from the isospin dependence of JV-JV cross section and the symmetry energy. 
Meanwhile, it is indicated that the squeeze-out flow depends strongly on the impact parameter and 
incident energy. Our study proposes that one can investigate the isospin-dependent reaction dynamics 
by studying the isospin effects on the azimuthal distribution and suggests that the isospin dependence 
of the squeeze-out flow could be as a probe of the isospin-dependent in-medium JV-JV cross section. 
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FIGURE CAPTIONS 

Fig. 1 The IDQMD model predicted Rn as a function of incident energy for reactions of ^^''Sn + 
^^''Sn and ^^^Ba + ^^'^Ba at impact parameter 6=4 fm. The lines are plotted to guide the eye. 

Fig. 2 The IDQMD model predicted Rn as a function of impact parameter for reactions of ^^"'Sn + 
^^''Sn and ^^*Ba + ^^''Ba at 350 MeV/nucleon. The lines are plotted to guide the eye. 

Fig. 3 The IDQMD model predicted normalized azimuthal distribution for ^^^Sn + ^^^Sn (solid 
circles) and ^^''Ba + ^^''Ba (open circles) at 350 MeV/nucleon and 6 = 6 fm for mid-rapidity 
nucleons by using different symmetry energy strength C and parametrizations of JV-JV cross 
sections: C = 32 MeV with experimental JV-JV cross section Ccxp (a), C = (no symmetry 
energy) with (Tcxp (b), and C = 32 MeV with Cugnon's JV-JV cross section crcug (c). Meanwhile, 
the results of Legendre polynomial fits according to Eq. (||) for ^^^Sn -I- ^^''Sn (solid line) and 
^■^''Ba + ^^"^Ba (dashed line) as well as the resulting ai and 02 are also included. 



TABLE CAPTIONS 

Table 1: The Rn values at different situations (see text) for i^'^Sn + ^'^'^Sn and ^'^^Ba + ^'^'^Ba at 
350 MeV/nucleon and 6 = 6 fm. 
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